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We investigate the thermodynamic phase transition taking place in the Blume-Capel model in 
presence of quenched disorder in three dimensions (3D). In particular, performing Exchange 
Montecarlo simulations, we study the behavior of the order parameters accross the first order 
phase transition and its related coexistence region. This transition is an Inverse Freezing. 



1. Introduction 



In the present paper, we will consider the Blume-Capel[ll-[1] model with quenched 
disorder (BC-random)[3-0]- ^ spin-glass model with bosonic spin— 1 variables 
{si = —1, 0, +1). BC-randorn^is one of the simplest spin glass models that displays 
an Inverse Transition (IT) [8]. By IT we mean a reversible transition occurring 
between phases whose entropic content is in the inverse order relation relatively 
to standard transitions. The case - already hypothesized by Tammann more than 
a century ago [9] - of "ordering in disorder" taking place in a crystal solid that 
liquefies on cooling, is generally termed inverse melting. The IT phenomenon also 
includes the transformation involving amorphous solid phases, as that of a liquid 
vitrifying upon heating, and the term inverse freezing (IF) is somewhat used in the 
literature: both phases are disordered but the fluid appears to be the one with least 
entropic content. IT has been experimentally observed in many materials: some ex- 
amples of inverse melting can be found in[lCH23l|. while IF takes place in 24l426l|. 
The reason for these counter intuitive phenomenon is that a phase usually at higher 
entropic content happens to exist in very peculiar patterns such that its entropy 
actually decreases below the entropy of the phase normally considered the most 
ordered one[27, 28]. 

The Mean Field (MF) solution of the BC-random model in the Full Replica Sym- 
metry Breaking (RSB) scheme predicts a phase diagram (fig.[T]) with a second 
order transition line (between Spin-Glass (SG) and Paramagnet (PM) phase) end- 
ing in a tricritical point, where a first order phase transition line starts, and from 
where a phase coexistence region departs. 

We stress that the transition is first order in the thermodynamic sense, with latent 
heat and is not related to the so-called random first order transition'^S^ occurring 
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in MF models for structural glasses. Furthermore, the first order transition is char 



acterized by the phenomenon of IF [30|,l3]|: the low temperature phase is PM, with 



a lower entropy than the SG phase, and the transition line develops a reentrance. 
In the present work we will study the behavior of the 3D BC-random model on a 
cubic lattice with nearest-neighbor quenched interaction. The nature of the phases 
that appear in the phase diagram, in particular across the IF First Order Phase 
Transition (FOPT), is studied through the shape of the order parameters distri- 
butions: this qualitative method allow us to understand, in a very simple way, the 
fundamental phenomenology that drives the IF scenario. Other results on the same 
model have been presented in 0, Is^] . 



2. Model and Observables 

The Hamiltonian of the BC-random is defined as follows 

nj[s] = -Y,JijSiSj + DY,sj (1) 

ihj) « 

where (ij) indicates ordered couples of nearest-neighbor sites, and Si = — 1,0, +1 
are spin— 1 variables lying on a cubic lattice of size N = with Periodic Boundary 
Condition. The external crystal field D plays the role of a chemical potential. 
Random couplings Jij are independent identically distributed as 

P{J^j) = - 1) + \KJ^, + !)• (2) 

We simulate two real replicas of the system and define the overlap, i.e. the order 
parameter usually characterizing the SG transition, as 

'^''^^Y.^^^) (3) 

i 

where (...) is the thermal average. If a thermodynamic first order phase transi- 
tion occurs, with latent heat, the most significant order parameter that drives the 
transition is the density p of magnetically active = 1) sites: 

P''' = ^Y^^si) . (4) 

i 

The apex J recalls us that the values of the parameters depend on the particular 
realization of disorder {{Jij})- All the information about the equilibrium properties 
of the system is in the knowledge of the following probability distribution functions 
(pdfs) 



P{q) = Pj{q) = 



P{p) = Pj{p) = 
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Figure 1. Left panel: the phase diagram of the BC model with quenched disorder for the fully-connected 
lattice with Gaussian distribuited random couplings (MF solution 5, 6]): the Second Order Phase Tran- 
sition ends in a tricritical point where an Inverse Freezing First Order Phase Transition takes place. The 
variance of P( Jij) in the MF model was oc l/z, z being the number of sites connected to each spin. In the 
3D model, where 2 = 6 (right panel), a bimodal distribution has been chosen with variance 1 rather than 
oc 1/6. 



We denote by ... the average over quenched disorder. The dependence on the 
random couphngs Jij is known to be self-averaging for the density probabihty 
distribution, but not for the overlap distributions Pj{q), whose average over the 
quenched disorder in the thermodynamic limit is different from the thermodynamic 
limit of a single realization of random couplings [s^]: 



P{q) = hm Pj{q) / lim > (q) . (6) 

We can introduce the notion of active and inactive site: when Sj = the site is 
inactive; otherwise, if = 1, it can interact with its neighbors (sj is an active 
site): as we will see in the next sections the IF indeed takes place between a SG 
of high density to an almost empty PM. The few active sites practically do not 
interact with each other but almost exclusively with inactive neighbors and this 
induces zero magnetization and overlap. The corresponding PM phase at the same 
D and high T has, instead, higher density and the paramagnetic behavior is brought 
about by the lack of both magnetic order (zero magnetization) and blocked spin 
configurations (zero overlap), cf. Sec 3. 



3. Phase Transitions and Order Parameters 



The equilibrium dynamics of the BC-random has been numerically simulated 
through Parallel Tempering (PT) technique: we have simulated in parallel the 
dynamic of the system at different values of T and D. For the PT in T, the swap 
probability of two copies at T and T + AT was: 

Pswap(A/3) = min [1, exp{A/3A?^}] . (7) 

While, for PT in D, two copies with D and D+AD were exchanged with probability 

Pswap(AZ)) = min [1, exp{/3AZ)Ap}] . (8) 

We will present data of 3D systems studied with PT in T at D = 0, and in D at 
T = 0.2, 0.3, 0.4, 0.5. At D = we simulated from 33 to 40 replicated copies Nt at 
linear size L = 6, 8, 10, 12 (number of disordered sample: Nj = 2000), for D = 
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Figure 2. The phase diagrams of the BC model with quenched disorder in three dimension. In the left 
panel is shown the T — D plane. The FOPT transition line develops a reentrance, it is an IT: decreasing 
the temperature the system leaves the SG phase and liquefies in a PM phase. In figure the spinodal lines 
(green) are reported and the blue line is the critical line of the FOPT. The arrow pointing to the region 
of the phase diagram where is located the tricritical point. In the right panel the T — p plane is shown. 



we simulated Nt G [20 : 33] at L = 16, 20 {Nj € [900 : 1500]) and Nt G [17 : 22] at 
L = 24 (TVj G [500 : 1000]). For the PT cycles in D, No S [21 : 37], parallel replicas 
at different D were simulated, for size L = 6, 8, 10, 12 and 15 {Nj = 1000). In the 
latter case, to resolve the coexistence region, varying AD were used, larger in the 
pure phases and progressively smaller approaching the transition. The number of 
Monte Carlo (MC) steps varies from 2^^ to 2^^ according to L and to the lowest 
values of T, D reached. 

Thermalization has been cross checked by looking at: (i) the symmetry of the 
overlap distributions Pn,j{q), (ii) the t-log behavior of the energy (when at least 
the last two points coincide), (iii) the lack of variation of each considered observable 
(e.g., P{q), Pip)) on logarithmic time- windows. 

In fig. ([2]) the phase diagram of the model is shown: both in T — D and in T — p 
plane. The second order phase transition has already been studied in jSj, details 
can be found in [32] . In short, by the study of the four-spins correlation function 



i i 

we can introduce a correlation lenght ^(T, L)/L that is scale invariant at the critical 
point: this property allows to define a size-dependet Tc{L). Finally, through Finite 
Size Scaling techniques, it is possible to calculate the critical temperature in the 
Thermodynamic Limit [33, 'ss']. 

The order parameter that drives the FOPT at finite L is the density distribution 
P{p)'- varying D,T, the system undergoes a transition with a discontinuous jump 
in p (and, thus, in q). The system is in the coexistence region if P{p) displays 
two peaks corresponding to the PM (low p) and SG (high p) phases. The first 
order transition line Dc{L,T), is the locus of points where the two phases are 
equiprobable, i.e., the areas of the two peaks are equal [sgI ]: 

dpPip) = / dpPip) (10) 

where po G [ppM '■ Psg] such that P{p) = (or minimal next to the tricritical 
point). In order to determine the transition point a method is to compare the areas 
under the distributions, cf. Eq (10). This is the point at which the configurations 
belonging to the SG phase and those belonging to the PM phase have the same 
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Figure 3. Density distribution Pi^{p), L = 15, across the coexistence region at T = 0.4 (left) and T = 0.5 
(right): two peaks develop at ppM and psG- As D increases the thermodynamically relevant phase (lovfest 
free energy) passes from SG to PM in a first order phase transition. The dominant phase correpsonds to 
the one with larger probability, i.e., larger intergral of the peak. As the peak at psG vanishes the system 
is in a purely PM phase. At T = 0.5 (right) a coninuum part between the SG and PM peaks appears: the 
distribution of the densities of two phases are overlapping. 



statistical weight: they yield identical contribution to the partition function of the 
single pure phase, and the free energies of the two coexisting phases are equal. We 
work at finite T moving D in that region and this method works quite well for 
T < 0.4 because the two peaks are clearly separated as soon as they appear, cf. 
fig. dS]). Determination of Dc is robust against reasonable changes of pQ. At T=0.5 
we have the problem that the distributions of the densities of the two phases 
are overlapping. In that case, seen the arbitrariness of choosing po, we actually 
determine the transition point as the D value at which the peaks have the same 
height. 

To have a better confidence with the results, two other methods have been used 
to determine the first order transition. These are not plagued by the problem of 
dealing with overlapping distributions since they do only rely on averages. The 
methods are: 

(i) Equal distance: at a given T, we plot D versus average p, we extrapolate the 
D{p) curves both from the PM and the SG phase {Dpm{p) and Dsg{p)) 
and we make a Maxwell-like construction determining a value of at 
which 



p{D,) = ]^{p^m{Dc) + Psg{Dc)) . (11) 

(ii) Equal area: equivalently (the equivalence is in the thermodynamic limit) 
one can determine -Dc as the value at which 

/ ppMiD)dD+ psGiD)dD= p{D)dD (12) 

JDsG JD, JDsg 



where I^sG find DpM are arbitrary, provided they pertain to the relative 
pure phases. The extrapolated, ppuiD) {psg{D)) is the inverse of the ex- 
trapolated curve Dpuip) {Dscip))- The curves are obtained by using all 
data, including those in the candidate coexistence region. 

We will show and compare in Sec. (3.2) the results obtained by these 
three methods. 
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Figure 4. PlC?) and Pl{p) for L = 16 and Nj = 1500 at D = 0.00. Across the transition P{q) (left 
panel) changes shape continuously from a Gaussian (in the high temperature phase) to, decreasing the 
temperature, a double peaked distribution with a continuum part between the peaks. In this region of the 
phase diagram the FOPT does not take place: P{p) (right panel) does not change shape. 



3.1. Second Order Phase Transition 

In fig. dH) (left panel) we report the distribution of overlap P{q) for a system of 
linear size L = 16 and D = 0: it changes shape accross the second order phase 
transition [s^] between the PM and SG phase from a Gaussian to a double peaked 
distribution 

Psciq) oc 5{q - qea) + S{q + qea) + f{q, L) , (13) 

where f{q,L) is a continuous function depending on the size. In the right panel 
of the figure we show the behavior of P{p) at fixed values of the crystal field 
and several temperatures: decreasing the temperature, deep in the SG phase, the 
average number of active sites is close to one 

hm >) = r dpP{p,T)p^l. (14) 

The distribution Pl{p) does not change shape across this transition. 



3.2. First Order Phase Transition 

Beyond the tricritical point, FOPT takes place and the system undergoes a dis- 
continuous transition between an "inactive" PM phase ((p) = ~ 0) and an 
"active" SG phase {pm 7^ 0). In the coexistence region, we can write the P{q) as 
sum of two contributes: 

P{q) = PsG{q) + PpM{q) . (15) 

For the PM contribution PpM^q) we have a Gaussian strongly peaked around q = 0. 
The Psciq) consists of a double peak (trivial) distribution with a continuum (non 
trivial) part between the two peaks, cf. Sec. 3.1. In fig. ([5]) we show the behavior of 
P{q) at different temperatures when the FOPT occurs. In the coexistence region, 
besides the double peak with a continuous part of the SG phase, a peak in g = 
appears due to the large density of empty sites. 

As shown in fig. ([2]), cf. also Sec. 2, FOPT occurs as an IF: the PM phase riches of 
inactive sites, below the T-range at which a pure SG phase is present, becomes the 
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Figure 5. P{q) across the FOPT at fixed T and different D values: we have a coexistence region between 
the phases SG and PM. The PM phase contributes to the P{q) distribution with a delta function in g = 0. 




low temperature (and less entropic) phase. This can be better observed in fig. ([6]) 
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1.98333(15) 
2.015(1) 
2.046(2) 
2.097(4) 


2.0243(95) 
2.043(5) 
2.092(5) 
2.143(4) 



Table 1. Critical values of the field D calculated through the equiprobabillty of the phases (eq. 9, blue line in 
the left panel of fig. |2]l), the equal distance (eq. 10) and the equal area (eq. 11) respectively. In the last two 
columns are reported the spinodal lines Dgp (green lines in the left panel of fog. (f2l'). 

where P{p) is represented at several values of D through the transition. In table ([T]) 
the critical values of the FOPT reported are obtained by the equal weight, equal 
area and equal distance methods. 

From the study of the shape of P{p) in the different phases (high-T PM phase, SG 
phase and low-T PM phase), the 3D BC-random model, allows us to interpretate 
the IF scenario in a very intuitive picture. Starting from a point of the phase 
diagram (D,T) in the high temperature phase (high-T and D> 0) the PM phase, 

as it can be seen in fig.Hl is dominated by the active sites with s? 7^ while {q"/"^) = 

(sp^sp^) = 0. Decreasing the temperature a Second Order Phase Transition takes 
place: P{p) does not change shape and pm{PM) = pm{SG). Decreasing further the 
temperature, the system undergoes to an Inverse-FOPT: P{p) becames a double- 
picked distribution. The coexistence of the phases occurs between a SG {pm > 0) 
phase and a PM {pm ~ 0) phase. Since the configuration space is dominated by 
inactives sites, the entropy of the PM phase becames smaller than the entropy of 
SG phase, leading to IF. 



4. Conclusions 

To conlcude, we analyzed the phenomenology of the BC-random in three dimen- 
sions studying the behavior of the order parameter distributions: P{q) and P{p)- 
In the case of the continuous transition, the lattice has an high density of active 
sites already the PM phaseEl 

lim f dpP(p,T,D)p r-. I (16) 
T^O,Z)=0 Jo 

and the density changes with continuity across the transition, whereas P{q) changes 
shape from a Gaussian to a double peaked distribution. When the discontinuous 
First Order Phase Transition takes place, P{p) becomes a double peaked distribu- 
tion due to the coexistence of PM and SG phases. The transition is an IF and 
the order parameter p jumps discontinuously between a low density and poor in- 
teracting PM phase (less entropic[§]) to an high density SG phase. Through the 
study of P{q) we have a clear evidence of the coexistence of a PM phase, riches 
in empty sites {pm 0) where the distribution has a single peak, and a SG phase 
with pm 7^ with a double peaked P{q) and a continuum part between the peaks. 

Finally we observe that the BC-random in three dimension is one of the few 
short-range disorder system undergoing a thermodynamic FOPT. We have found 
in literature only one system that displays a FOPT: the 4— Potts glass studied by 



Fernandez et al. [381]. In that case, though, randomness tends to strongly smoothen 



the transition into a second order one and the finite size effects are very strong in 



^In the limit limy.^o,l5~+-oo Jq^ dpP(p,T, D)p = 1 and we obtain the Edward-Anderson model [40ll . 
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determing the tricritical point. Our study, thus, confirms in a clear way the claim 
of the existence of a FOPT in presence of quenched disorder thanks to almost neg- 
ligible finite size effects. 

Moreover, we notice that the FOPT of the BC-random is driven only by exter- 
nal thermodynamics parameters, e.g., the temperature and the chemical potential. 
Even though, from the point of view of the numerical simulation changing the 



pressure, bond dilution [39(| or even the relative probabilities of having ferro- or 
antiferro- magnetic interactions [1^ is technically equivalent, the latter are com- 
plicated to control in a real experiment and require the preparation of several 
samples with different microscopic properties. Eventually, we are not privy to any 
three dimension short-range system with quenched disorder undergoing an IF. 
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